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We present a dynamical analysis of the inflationary flow equations. Our technique uses the Ifubble 
'jerk' parameter ^ = {mpY/16iT^)H~^ {dH/dcj>d^H/d(f)'^) (where H is the Hubble parameter and <jf> 
the inflaton) as a discriminant of stability of fixed points. The results of the analysis are used to 
explain qualitatively the distribution of various observable parameters (e.g. the tensor-scalar ratio, 
r, and scalar spectral index, Us) seen in numerical solutions of the flow equations using stochastic 
initial conditions. Finally, we give a physical interpretation of the flow in phase-space in terms of 
slow-roll motion of the inflaton. 
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1. INTRODUCTION 

Since the first discussions of inflationary cosmology 
more than twenty years ago, a huge number of spe- 
cific models have been proposed. Most of these are 
'phenomcnological', though increasing effort has been in- 
vested towards constructing physically motivated models 
such as brane inflation (for reviews, see [J, Q)- Given a 
plethora of models, one can ask whether aspects of in- 
flation can be studied without recourse to particular sce- 
narios. 

Hoffman and Turner [sj first addressed this by intro- 
ducing flow equations that govern evolution of observ- 
ables in any single-field model of inflation. These flow 
equations were later generalised to arbitrary order by 
Kinney In his notation, the inflationary flow equa- 
tions are : 



de 
dN 
da 
dN 
d i. 
dN ^" 
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Here the derivative with respect to the number of e- folds, 
N, runs in the opposite direction to time. The flow equa- 
tions represent an infinite dimensional dynamical system 
of the following parameters given in terms of the Hubble 
parameter H{(j)): 
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where mpi is the Planck mass and primes denote deriva- 
tives with respect to the inflaton cj). The flow parameters 
{e, (T,^ Aijjare related to the conventional slow-roll pa- 
rameters (a| 



'"PI 



(3) 



etc. defined in terms of derivatives of the inflaton poten- 
tial V{4)). The validity of the slow-roll approximation is 
more easily expressed in terms of these parameters than 
those of equation The relationship between the two 
sets of parameters is discussed in §3. 

Solving the flow equations gives e{<j)) or, equivalently, 
H{(j>). These in turn fix the shape of potential V{(l)) via 
the Hamilton-Jacobi equation 
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which is equivalent to the Friedmann 'acceleration' equa- 
tion 



(5) 



From ([5]), it is clear that inflation is sustained as long as 
e < 1. Recently, Liddle showed that an analytic solu- 
tion to the flow equations, closed at some order, yields a 
potential V{(f>) expanded as a rational polynomial. 

Using the flow equations, a large number of numeri- 
cally viable inflationary models can be readily produced 
stochastically. However, a mathematical analysis of the 
flow equations - an infinite hierarchy of differential equa- 
tions - appears daunting. While there have been discus- 
sions of 'attractor' solutions l7| in the past, in this 
paper we present a more rigorous analysis of the fiow 
equations. We begin by exploring the dynamics in the 
phase-space of flow parameters {e,a,^ Xh} and classify- 
ing all trajectories. In particular, we are able to explain 
the complicated motions in the (cr, e) plane as observed 
in [1, 0] ■ We then analyse how these motions translate 
to structures in the space of observables. Finally, we dis- 
cuss the physical implications of our analysis when the 
Hamilton-Jacobi equation is incorporated. 
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2. PHASE-SPACE DYNAMICS 
2.1. The nature of fixed-points 

Two sets of fixed points of the flow equations are evi- 
dent from inspection: 

• e = 0, (J — constant- so called 'r — 0' fixed points, 
from the fact that the tensor-scalar ratio r cx e 



vanishes along that line. 



fixed point, 
power-law inflation from 
V{(j)) oc exp ( — \/87r(/)/mpi). 



• e > 0, cr = -2e 

'power-law 



^Xh)- so called 
being a feature of 
a potential such as 



As a starting point, let us examine the nature of these 
fixed points in the simplest two-dimensional projected 
phase-space. 



2.1.1. Simplest [a, e) dynamics 

Consider the two lowest order equations in the system 
([2]) . The dynamics in the (ct, e) phase-space clearly de- 
pend on the value of ^Aff, the Hubble 'jerk' parameter. 
Following the notation of we shall rewrite it as 
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Thus we are examining the dynamics of the two- 
dimensional system 



dN 
da 

dN 



= e(fT + 2e) , 

= -e(5CT + 12e) + 2^ 



(7) 
(8) 



The fixed points in the (cr, e) system are easy to analyse 
if ^ is treated, for now, as a constant. By linearising 
d?]) and (HI) , one sees that the sign of ^ discriminates the 
stability of fixed points in the (cr, e) plane. The results of 
a standard fixed point analysis is shown in Table I. 

We draw attention to the 'nuUclines' where the flow 
becomes vertical or horizontal: 



Horizontal nuUcline 
Vertical nuUcline 
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(9) 
(10) 



These curves provide useful demarcations of where flows 
change direction, and we include them in our phase por- 
traits. 



2.1.2. Methodology 

Inflation can be investigated 'stochastically' by first 
selecting a random point in the flow parameter space as 



an initial configuration of a particular universe. Because 
this is not the same as assigning randomness to initial 
conditions for inflation, we adhere to the phrase 'initial 
conflguration' to keep this distinction clear. 

An initial conflguration is then evolved forward in time 
(backward in e-fold) until inflation ends in one of the 
following plausible ways: 

1. By achieving e = 1. When this happens, we say 
for convenience that the 'slow-roll' condition has 
been violated (keeping in mind that inflation ends 
because equation ^ changes sign rather than be- 
cause the slow- roll parameters ey and \riv\ of equa- 
tion ([3]) exceed unity). Observables are then cal- 
culated a specific number of e-folds (which we take 
to be 55, unless otherwise stated) before the end 
of inflation by integrating backward in time. This 
choice is in accordance with analyses of references 
[1, which indicate that for plausible assumptions 
concerning reheating etc, the observable perturba- 
tions were generated 50 to 60 e-folds before the end 
of inflation. Variation of this number is explored in 
§2. 

2. By an abrupt termination, perhaps from interven- 
tion of an auxiliary field in hybrid inflation or 
when open strings become tachyonic in brane in- 
flation Because these scenarios accommodate 
large number of e-folds during inflation, one identi- 
fles them with an asymptotic behaviour of a trajec- 
tory. Observables are calculated along the asymp- 
tote. We return to the practical issues of calculat- 
ing observables later. 

Because we evolve flows backward in e-fold, the no- 
tions of stability and attraction found in Table I become 
reversed. This is summarised in Table H, from which we 
conclude that given an initial configuration in {e, Xh} 
space, as long as the function: 



sgn(^) 



remains constant, then the dynamics in the {e,a,^ Xh} 
phase-space are completely characterized. 



2.1.3. Dynamies of ^ 

Having assumed that ^ was constant, we can now in- 
corporate aspects of the dynamics of ^ using the function 
sgn(^). If sgn(^) changes unpredictably, then the classi- 
fication in Table II would be useless. We now show that 
this is not the case. 

Starting from equation one finds that the evolution 
of ^ has no explicit dependence on e 




d^ 
dN 
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TABLE I: Stability of fixed points in the (tr, e) plane with respect to A^, 
the number of e-folds before inflation ends. Vertical and horizontal nuUclines 
are included in the phase portraints. In the case where ^ = 0, flows coming 
from the right become tangential to the luie a = ~3e, which is shown. 





Attractors 


Comments 


Sketch 




{e = l,c7 <-2} 


- Ail trajectories eventually vio- 
late sbw-roll. 

- No asj mpcote 


"^^^^^^ 


£ 

T' = o 




{f = l,cr < -2} 
{e = 0,0- > 0} 


- Reached by flows starting on the 
left of the critial trajectory T 
passing through (—2, 1) 

- flow starting on the right 
of T approaches this point 






^ <0 


{e= l,cr < -2} 
{e = 0,cr ^ .do} 


- same as the case <^ = 

- flow starting on the right 
of r approaches this line 


' ^ ^ =2 









TABLE II: Stability in the (a, e) plane with respect to time with condition e < 1 imposed. 
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The solution of equation ((TT|) is 



2.1.4- Numerical method 



N 



a(n)dn 



No 



^(7V)=e(A^o)exp(i 

+ / ^AH(»7i)exp( / a{n)dn\dm . (12) 

J No ^ '2 Jn ' 

Now we use equation to deduce the following points: 

• If a trajectory does not violate slow-roll (i.e. in- 
flation ends abruptly), then from Table II, there 
exists some e-fold A'o along such trajectory such 
that C(iVo) < and cr(7Vo) > for all iV < iVo 
(since trajectories follow dN < 0). Assume that 
typically \^Xh\ ^ £. and using the fact that a{N) 
is unbounded for asymptotic flow, yields 



hm C = 0_ 

N^-oc 



(13) 



In other words, if ^ starts out negative, then it 
is exponentially damped towards zero from below. 
Slight deviations may be expected if "^Xh is large, 
but this is atypical. 

• If inflation ends at e = 1 , a similar argument (tak- 
ing care of the different integration direction) shows 
that 



lim C = 0_ if ^{No) < . 

N^-\-oc 



(14) 



If ^(iVo) > 0, we also have exponential decay to- 
wards zero, but only for 55 e-folds. Taking Nq = 
to be the end of inflation, we have : 



/I 

e(55)=e(0)exp(- a{n)dr, 

/.55 / 1 /■"' \ 

+ / '^'^^("i) exp ( — - / a{n)dn\dm , (15) 
^0 V 2 J55 / 



at 55 e-folds prior to the end of inflation. Un- 
like the asymptotic cases, ^ cannot have decayed 
by very much after only 55 e-folds, and indeed the 
magnitude of ^(55) can be rather significant espe- 
cially when, say, > 0. Typically one finds 
^(55) - 0(10-2). 

The conclusion is that regardless of its initial value, ^ 
always exponentially decays towards zero before observ- 
ables are calculated, thus typically preserving sgn(^). In 
particular, if inflation ends with ^ > 0, then the value 
of ^ when observable perturbations were produced may 
be significantly non-zero and positive. We return to this 
point when we study the structures in the space of ob- 
servables. 



We ran a computer code (previously used in [l^]) that 
follows the stochastic prescription to see if we can un- 
derstand actual phase-space trajectories. The program 
selects initial configurations at random from uniform dis- 
tributions within the following ranges : 

eo e [0,0.8] , 
ao G [-0.5,0.5] , 

^0 e [-0.05,0.05] , (16) 
'^XhIo e [-0.025 X 5-^+^ 0.025 X 5-^+3] , (3 < £ < 10) 
"A^flo = , 

where the truncation is imposed to close the hierarchy at 
£ = 10. For most applications of the flow equations, the 
exact widths of distributions in (jl6p are unimportant as 
long as the hierarchy is convergent. Figure 1 shows (tr, e) 
trajectories from nine initial configurations. They are all 
consistent with the linearised analysis in Table II. 

By producing 125,000 realisations, we observed the fol- 
lowing distribution among three classes: 

(A) Asymptotic flow towards the half-line e = 0, cr > 
(for « 90% of all trajectories). 

(B) Slow-roll violation with observable perturbations 
produced at the power-law line a = — 2e (for « 8% 
of trajectories). 

(C) Slow-roll violation with observable perturbations 
produced at e = 0,cr < (for 2% of trajecto- 
ries). 

This is in agreement with observations in the numerical 
investigations of references [1, [l^, [IH . We have a simple 
explanation for this distribution. 

Figure 2 shows the sketches of Table II with the initial 
window (jl6p in the (cr, e) plane superimposed, and regions 
colour-coded according to the eventual fate (A, B or C 
above) of points starting within them. One deduces that 
the ratio of these colour-coded areas determines the final 
distribution in a given stochastic experiment. Formally, 
for the uniform distributions of equation (16) we have : 



Percentage of points in class X 
area for outcome X 

area of (an, en) window 
sgn(4) V u, u/ 



(17) 



where X = A, B or C described above. The sum can 
be evaluated for ^ = 0,±(|^o|) ~ 0,±max(^)/2. Apply- 
ing (flT)) to our case where ±(|Co|) — ±0.025, one easily 
deduces from Figure 2, for example, that the dominant 
yellow area in the three panels corresponds to the ma- 
jority of points belonging to class (A), whereas the tiny 
cyan areas mean that models belonging to class (C) are 
very rare indeed. 

Figure 2 also shows the effect of changing the size 
and location of window ((To,eo) in (|16p . For instance. 
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FIG. 1: Sample of nine (cr, e) trajectories. Circles denote initial random points. Crosses mark where the observable perturbations 
are produced. In each panel the blue line is the horizontal nuUcline ® Models A,B,C show the most generic behaviour for 
all values of sgn(^) (see Table II). The remaining models all violate slow-roll and observables are calculated by backward 
integration. According to our classification, models D,E and F are clearly consistent with ^ > 0, whereas ^ < for G, H and I. 



^=0.025 




FIG. 2: Illustrations of critical areas in the {a, e) plane. A random point that starts in a yellow region ends up around 
e — 0,a > (class A in text). Pink and cyan mark potentially slow-roll violating points whose observables are produced at 
a = — 2e (class B) and e = 0,ct < (class C) respectively. Boxes indicate the initial (cro,eo) window in equation (|16p . The 
coloured areas translate to respective abundance of resulting points from a stochastic experiment as summarised in equation 
(|17[1 . However, this breakdown is sensitive to the parameter A'critical (^^® number of e- folds before infiation is terminated 
abruptly, see text for a more precise definition). For instance, very large A'^^.j^j^g^j renders the top panel entirely pink. 
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if max(eo) < 0.8, the percentage of points in case (A) 
will drop proportionally. The fraction of trajectories in 
each of the classes (A), (B) and (C) depends primarily on 
the initial distributions of the parameters a and e, rather 
than the initial distributions of higher-order parameters. 

Finally, we comment on a particular aspect of numer- 
ical solutions of the flow equations. Consider the case 
when ^0 > 0. From Table II, one expects all trajecto- 
ries to violate slow-roll. In particular, a trajectory may 
appear to approach the line e = 0, cr > 0, but even- 
tually it inevitably intersects the vertical nuUclinc (fTU)) 
and swiftly reverses direction, thereafter violating slow- 
roll. In practice, however, one identifies a flow to be 
asymptotic (and obeying slow-roll) if it supports infla- 
tion until a prescribed critical number of e-folds is ful- 
filled. In our code, this number is set to A'criticai = 200. 
Therefore, trajectories that could eventually violate slow- 
roll for N > A^criticai SLTC uot identified as such in prac- 
tice. This explains why in Figure 2 there are points that 
will not violate slow-roll even though ^ > 0. Of course, 
to select only truly asymptotic flows, one may increase 
A^criticai to bc, Say, 1000 or more. However, this is at the 
expense of proportionally increased computation time. 

To summarise, the numerical implementation of the 
flow equations described here agrees with theoretical 
analysis of the flow dynamics, based on the parameter 
^, described in the previous Section. There is some de- 
pendence of the numerical solutions on input parameters 
such as the e-fold cut-off A'criticai and choice of initial 
(cr, e) window (with weak sensitivity to distributions of 
higher order flow variables). These dependences can be 
understood theoretically. In particular, we interpreted 
the distribution of numerical outcomes as ratio of areas 
in the theoretical phase-space. 



3. STRUCTURES IN THE SPACE OF 
OBSERVABLES 

At about 55 e-folds before the end of inflation, ob- 
servable perturbations were produced, setting the values 
of key inflationary observables- for example, the scalar 
spectral index Us , tensor-scalar ratio r and running of 
spectral index -jj^- To second order, these are related 
to the inflationary flow variables by [13j : 



r ~ 16e[l - C(cr -H 2e)] , 

T,, ~ 1+0" 



(18) 

(5 „ 3C)e2 _ 1 (3 _ 5C)ae + 1(3 - C)^ , 

(19) 
(20) 



d\nk 1-e dN ' 

where C = 0.0814514. Note that define r to be 16 

times smaller than in equation (|18p . 

We can now make the following classification: 

• If inflation ends by violation of slow-roll conditions, 
then according to our earlier analysis, observable 
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FIG. 3: Models in the {ns,r) plane generated from the flow 
equations. The line r = is densely populated. The swathe 
of points traces the curve given by (|21|l (shown by the blue 
line) . The distribution of points along this line is discussed in 
text. 



perturbations were produced either at (i) the line 
e — — 2(T, or (ii) e — 0,a < 0. In the flrst case, we 
have an exact power-law inflation, in which case the 
relation 



1 



2r 



16 



holds f3|. For case (ii), one observes 







< 1 



(21) 



(22) 



• If inflation does not violate slow-roll but ends 
abruptly, then we must have asymptotic flow along 
e — 0,a > 0. This yields a 'blue' spectrum : 







> 1 



(23) 



Equations ((2T|) . ((22|) and ((23|) deflne three prominent 
structures in the (n^ , r) plane where we expect points to 
cluster strongly. From the analysis of the previous sec- 
tion, we also know the relative density of these structures: 
Most points cluster at (class A) , some at the power- 
law curve pT|) (class B), and very few at ([22|) (class C). 
This behaviour is seen clearly in Figure 3, which shows a 
scatter plot from 125, 000 numerical solutions of the flow 
equations using the code and initial conditions described 
in the previous section. 

In terms of the parameter ^, we deduce that the swathe 
corresponding to equation (|2ip is made up almost en- 
tirely of models in which ^(55) > 0. The r = clus- 
ter should theoretically be made up of only those mod- 
els with ^(55) < 0, although in practice, the choice of 
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itical 



means that there are some points with ^(55) > parameters : 



on this line as well. 

The swathe is not evenly populated along the curve 
(PT|) . We have a simple explanation for this non- 
uniformity. Starting at the end of inflation with ^(0) > 
and e(0) = 1, the trajectory is traced 55 e-folds backward 
in time, at which point the trajectory must be near the 
line (7 = — 2e in order to form the swathe in the (ns,r) 
plane. From Table I, the e coordinate of the power-law 
fixed point is 



From (US]), to first order, 



16e(55) > 16e* = 16a/^(55) 



so that 



16^/£,{55) < r < 16 



(24) 



where ^(55) is given by equation (|15|) . Only points within 
these bounds are expected to lie close to the line defined 
by equation (^J). Intuitively, we should expect to see 
some sharp features in the swathe because trajectories 
that give rise to the swathe are suddenly terminated 55 e- 
fold by backward integration in order to calculate r. For 
example, in Figure 3, the lower bound in (|24p is found to 
be 1.7 on average. For r < 1.7, the points show a large 
scatter in the {rig , r) plane and no longer cluster tightly 
around the relation defined by equation l[2T|). 

Finally, we comment on the sensitivity of the swathe to 
the e-fold limit iVobs = 55 when observable perturbations 
are generated. As discussed earlier, the investigations 
of [E favour a range of A'^obs between 50 to 60, but 
caution that some models with prolonged reheating may 
give A^obs as high as 100. We investigate how such an 
increase affects the structures in the {ug^r) plane. 

An increase in A"obs only affects slow-roll violating 
models. In the case where observable perturbations are 
produced along the line e = 0, (t < 0, extended back- 
ward integration gives rise to a larger spread along the 
line ris < 1 in the observable plane. On the other 
hand, for those trajectories approaching the power-law 
line a — —2e, an increase in A'obs means that ^ exponen- 
tially decays towards zero for a longer time, concentrating 
the points to lower values of r. Furthermore, as A"obs is 
increased, the analogous lower bound to ([24]) extends to 
lower values of r. Thus for large values of N^hs, points 
cluster tightly around equation (|2ip to values of r close to 
zero. These effects are evident in Figure 4, which shows 
analogous plots to Figure 3, but for numerical simula- 
tions with A"obs is the range 100 to 500. 



4. PHYSICAL INTERPRETATIONS 

In the slow-roll approximation, the potential shape is 
parametrized by its slope and curvature via the slow-roll 
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(25) 



We wish to see precisely how the flow parameters 
{e,a/ Xh} relate to this description of the potential 
shape. 

As stated in the introduction, the flow parameters are 
related to ¥{(/)) via the Hamilton- Jacobi equation ([4]). 
Dividing it by gives a useful form : 



(26) 



By differentiating we can derive two further relations: 

(27) 



772-3(e + r;)+e = - 



(28) 



where the sign of ^/e matches that of H' . Equations ([26]) 
to relate the amplitude, gradient, and curvature (in 
that order) of the potential V to the flow parameters. 

Now consider the flrst slow-roll condition ey ^ 1. This 
condition states that the slope of the potential should 
be so small that (i) potential energy dominates kinetic 
energy of the field {V (fP') and (ii) rolling approaches 
'terminal' speed so that « 0. By dividing ([?7)l by 
we obtain these requirements exactly in terms of flow 
parameters: 



^-3 
e- 3 



< 1 



(29) 



One way to maintain inflation is to ensure that the po- 
tential becomes increasingly flat, so that the inflaton con- 
tinues approaching the minimum asymptotically. From 
this can be achieved with a small e regardless of 
how large is. In the (cr, e) phase-space, this corre- 
sponds to those 'hybrid' type trajectories approaching 
e = regardless of large a. 

Next, the second slow- roll condition |77y| ^ 1 can be 
obtained by differentiating ey, and so it places a restric- 
tion on the potential curvature. By dividing by (|26|) . 
we can recast this condition as 



3(e + r;) - 772 - e 



< 1 



(30) 



(Note the role of parameter ^ as a constraint on the cur- 
vature.) In particular, we see that a universe with suffi- 
ciently large negative 77 will violate ([50]) even though the 
first slow-roll condition ((29)) is satisfied. This shows why 
potentials with sufficiently large negative curvature (such 
as those associated with spontaneous symmetry breaking 
or new inflation [l3,[lB|) are doomed to violate slow-roll. 
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FIG. 4: Numerical simulations similar to those shown in Figure 3, but for different values of A^obs- A^obs increases, the 
swathe of points around equation lpT)l moves to lower values of r as described in the text. 



As the inflaton gains speed along a concave-down poten- 
tial, its 'image' in the phase-space is a flow towards the 
line {e = 1, cr < —2} as shown in Table II. 

In summary, one can understand the correspondence 
between inflaton motion along a potential and the flow in 
{e, cr,^ A//} phase-space simply by using the flow param- 
eters to characterize the gradient and curvature of the 
potential. 



5. CONCLUSIONS 

The inflationary flow equations have been used by a 
large number of authors to draw inferences about ob- 
servable parameters. When solving the flow equations, 
one finds favoured regions of observable parameters due 
to dynamical attractors. The purpose of this paper has 
been to examine the dynamics of the flow equations an- 
alytically. In particular, we have shown that the Hubble 
'jerk' parameter effectively discriminates the stability 
of fixed points in the flow parameter phase-space, while 
the dynamics of ^ itself can be described simply in terms 
of exponential decay. Using this technique, we were able 
to explain complicated motions in the {e, cr,^ Xh} phase- 



space analytically. The distribution of points and stabil- 
ity of structures in the plane of observables can also be 
understood in terms of dynamics of ^. Comparing our an- 
alytic work with numerical solutions of the flow equation, 
we found that they were in good agreement. Finally, we 
have given exact relations showing how the trajectories 
in the phase-space of flow variables correspond physically 
to motion of the inflaton along a potential. 

The main features in the plane of observables are gen- 
erally insensitive to the initial configuration (i.e. choice 
of values ^A/f|p) - although there is a dependence on the 
parameter iVcriticai which can be understood dynamically. 
Thus we believe that we understand the main features of 
the solutions of the flow equations. Nevertheless, to inter- 
pret diagrams such as figures 3 and 4 in any statistical 
sense requires an understanding of the thorny issue of 
measures and how initial conditions for inflation might 
translate into distributions of initial configurations. We 
have not attempted to address these problems in this pa- 
per. For some speculations, see references [igI [iTj. 
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